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Abstract 

A thermodynamic argument is proposed in order to discuss the most appropriate form of 
the local energy balance equation within the Oberbeck-Boussinesq approximation. The 
study is devoted to establish the correct thermodynamic property to be used in order 
to express the relationship between the change of internal energy and the temperature 
change. It is noted that, if the fluid is a perfect gas, this property must be identified 
with the specific heat at constant volume. If the fluid is a liquid, a definitely reliable 
approximation identifies this thermodynamic property with the specific heat at constant 
pressure. No explicit pressure work term must be present in the energy balance. The 
reasoning is extended to the case of fluid saturated porous media. 
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1 Introduction 

The Oberbeck-Boussinesq approximation, so named after the pioneering works by Oberbeck 
PQ and Boussinesq [2], is the basis of most of the contemporary studies on natural or mixed 
convection flows. Very interesting historical surveys on the origins of this approximation are 
available in the recent papers by Zeytounian [3] and Bois [JJ. 

Although the nature of this approximation is very clear and unambiguous with reference to 
the mass and momentum balance equations for the fluid, the formulation of the approximated 
energy balance equation is not so definite and univocal. The questions concerning the energy 
balance equations are the following: 

A) Which is the specific heat involved in the energy balance? 

B) Is there a pressure work term in the energy balance, proportional to the convective 
derivative of the pressure field? 

The textbooks on fluid dynamics and heat transfer generally give clear answers to these 
questions. The problem, as it will be discussed in Section [3l is that the answers are different. 

In a recent technical note [5] , the present author carried out a first analysis of the existing 
formulations of the local energy balance adopted in the Oberbeck-Boussinesq approximation 
of buoyant flows. 

The purpose of this short paper is to extend the analysis performed in Ref. [5] in order to 
point out the manifold nature of the energy balance formulations, within the framework of 
the Oberbeck-Boussinesq approximation, available in the literature. Then, a thermodynamic 
argument is proposed in order to give answers to questions A) and B). In particular, it will be 
concluded that the answer to question A) depends on the fluid being a liquid or a gas. For a 
perfect gas the answer to question A) is definitely: "the specific heat at constant volume c^" . 
For a liquid the answer to question A) is less definite, but sufficiently reliable: "the specific 
heat at constant pressure c p ". The answer to question B) is: "no pressure work term appears 
in the energy balance" . In a final section, the analysis of the energy balance is applied to 
the topic of buoyant flows in fluid saturated porous media. 
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2 Mass and momentum balance 

As is well known, the Oberbeck-Boussinesq approximation implies that the local mass and 
momentum balance equations be written as 

DUi 1 dp e 2 

-fTT = ^ {T -Tq) [5gi + vV Uil (2 

Dt po oxi 

where the summation over repeated indices is assumed. In Eq. (J2]), D/Di is the substantial or 
convective derivative. In Eqs. ([I]) and ([2]), u j is the velocity field, Xi is the position vector, t is 
the time, T is the temperature, gi is the gravitational acceleration, v is the kinematic viscosity, 
(3 is the isobaric coefficient of thermal expansion, po and To are the reference density and the 
reference temperature respectively, while D/Di is the substantial or convective derivative. In 
Eqs. ((T]) and ©, the properties v and (3 are also referred to the temperature To. In Eqs. (pQ) 
and ([2]) , the properties v and (5 are referred to the temperature To . The implicit assumptions 
behind Eqs. ([T]) and ([2]) are that: p e = p — po gi x^ is the difference between the pressure 
p and the hydrostatic pressure, and that one considers the density as coincident with the 
reference value po except for the gravitational body force term. For that term, the density 
p is assumed to be a function of the temperature only, thus considering the dependence on 
the pressure as negligible. The linear equation of state 

p(T) = po[l-/3(T-T )] (3) 

is implicitly invoked in Eq. $2). In Eq. ([3]), the dependence on T is assumed to be sufficiently 
weak to be approximated linearly in the surroundings of the reference value To. 

3 Energy balance 

For the energy balance, the formulation of the Oberbeck-Boussinesq approximation is not 
so definite in the literature. In fact, one may have Chandrasekhar's [6] and White's [7] 
formulation 

DT 

p c v — = kV 2 T + 2p® lj ® lj , (4) 
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where k is the thermal conductivity, p is the dynamic viscosity and T>ij is the strain tensor 

The source term in Eq. ([5]), 2 pDijDij, is the thermal power generated by the viscous 
dissipation. 

One may have the enthalpy formulation [HJ [§] 

Poc p — = k\7 2 T + 2pV lj V ij + (3T^ 1 (6) 

where the last term on the right hand side is an additional source term: the pressure work 
acting on the fluid element. 

Finally, one may have Landau-Lifshitz's [10j . Bejan's and Kundu-Cohen's [T2] formu- 
lation 

DT 

p c p — = k V 2 T + 2 MJDij % • (7) 

In order to decide on the most convenient expression of the energy balance, let us write the 
general non-approximated form of this balance [6JQ3], i.e. the local version of the First Law 
of thermodynamics 

p m = -^ + a ^ (8) 

where e is the internal energy per unit mass, cr j,- is the fluid stress tensor and qi = — k dT/dxi 
is the heat flux density. The meaning of the terms on the right hand side of Eq. dSJ) is 
straightforward: the first term is the incoming heat flux contribution to the energy change, 
while the second term is the mechanical work input due to the forces acting on the boundary 
of the fluid element. The latter term depends only on mechanical quantities, i.e. on the 
velocity and pressure fields within the fluid domain. On the other hand, the evaluation of 
the energy change on the left hand side of Eq. ([HJ implies the use of thermodynamics. 

Thermodynamics ensures that e = e (T, p) for every single-phase or two-phase stable equi- 
librium states. In the special case of a perfect gas, it is well known that e = e (T) [JJ], so 
that 

de = c v dT. (9) 
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In the case of either a liquid or a real gas, one must rely on the Oberbeck-Boussinesq ap- 
proximation by assuming that an approximate equation of state p = p (T) can be applied. 
This implies that the pressure of the fluid does not change appreciably. Since p = p (T) and 
since the pair (T, p) yields a unique stable equilibrium state, then one concludes that all the 
thermodynamic properties may be considered as functions of T. This conclusion holds for 
the internal energy per unit mass, so that a relationship 

de = cdT (10) 

can be established. The thermodynamic coefficient c, in general, does not coincide either 
with c v or with c p . In fact, c is the total derivative of the function e = e(T,p(T)) with 
respect to T, and not the partial derivative of e = e (T, p) with p kept constant. As is well 
known, the latter is the correct thermodynamic definition of c v [13]. The equation of state 
p = p (T) is one regarding a set of stable equilibrium states of the fluid with approximately 
the same pressure. Then, one has 



de 
df 



(11) 

v 



Eq. (jlip is not the definition of the specific heat at constant pressure c p . As is well known 
the latter is defined as 

<*),• (12> 

where h = e + p/p is the enthalpy per unit mass. Then, one can easily write the following 
relationship: 

c = c p , (13) 

P 

where the definition of the coefficient of isobaric expansion 

has been used. Then, c is smaller than c p and differs from c v , except for the limiting case 
of a perfect gas. Indeed, in the latter case, one can easily show that the equation of state 
of the perfect gas and Eq. (fT3l) ensure that c — c v , so that Eqs. ([9]) and (fT0|) are perfectly 
consistent. 
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T[°C] 



Figure 1: Comparison of c with c p and c v for water at atmospheric pressure 



8. xlO" 6 




Figure 2: Comparison of c with c„ for water at atmospheric pressure 



With reference to a Newtonian fluid and on account of Eq. (TjQ), one may express o"y 3?« 
2 fiDij Then, by employing Eq. (fTUI) . by replacing p with po an d by assuming k 
constant, Eq. ([8]) can be approximated as 
DT 



(15) 

In Eq. (|15p . the value of the thermodynamic coefficient c, given by Eq. (|13p . is referred to the 
temperature To like any other fluid property involved in the Oberbeck-Boussinesq system of 
governing equations. Eq. (|15p coincides with Eq. (jlj) only for the limiting case of a perfect 
gas and differs from Eq. (J7|) as c < c p . 

One can question the extent to which c p and c v differ from c in the case of liquids. For water 
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Table 1: Specific heat ratios for some organic liquids at 25° C and atmospheric pressure 





c)/ Cp 


(C CyJ / ' C v 


[Cp Cv) 1 Cp 


Acetone (C 3 H 6 0) 


0.0087% 


40% 


29% 


Benzene (C^Hq) 


0.0075% 


36% 


26% 


Ethanol (C 2 H 6 0) 


0.0074% 


37% 


27% 


Methanol (CH 4 0) 


0.0076% 


37% 


27% 



at atmospheric pressure, a precise evaluation of the discrepancies (c p — c)/c p , (c — c v )/c v and 
(c p — c v )/c p can be done by means of the data reported in Appendix C of Bejan's textbook 
[llj. In fact, on account of Eq. (|13j) . one has 

Cp-c pj3 c-c v p(c p -c v )-p(3 p(c p -c v )-pP 

= , = = — — r. (16) 

Cp p Cp C v p C v p Cp p yCp C v ) 

Figure [1] evidently suggests that the discrepancy (c — c v )/c v is rather close to (c p — c v )/c p 
and that (c p — c)/c p is much smaller. Figure [2] shows that (c p — c)/c p is, on average, of order 
10~ 5 . As a consequence, for water at atmospheric pressure, the approximation 

c ^ Cp (17) 

is a definitely reliable one. In general, it is not easy to find data for the specific heat at 
constant volume of a liquid. Usually, thermodynamic tables report the values of c p , while c v 
is evaluated from the Mayer relationship [J3] 

pK T 

where kt is the coefficient of isothermal compressibility, 



(18) 



On account of Eqs. (|16p and (]18p and of the data reported in Refs. |14|. I15j . one may evaluate 
the discrepancies (c p — c)/c p , (c — c v )/c v and (c p — c v )/c p for some organic liquids at 25° C and 
atmospheric pressure. These data are reported in Table HJ This table suggests again that 
Eq. (|17p is definitely reliable with an error smaller than 0.009%, while the approximation 
c = c v would lead to an error higher than 35%. 
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4 Possible pitfalls 

In the preceding section, a thermodynamic strategy has been established to determine the 
most appropriate formulation of the local energy balance equation within the Oberbeck- 
Boussinesq approximation. The basis of this approach is twofold. 

• The convective derivative of the internal energy per unit mass, De/Di, is evaluated by 
considering the thermodynamic process undergone by the fluid element. This process 
can be reliably modeled as an isobaric process. 

• The mechanical work input er^ Dij is evaluated according to the stress-strain relation- 
ship for a Newtonian fluid, as well as to the constraint dui/dxi = satisfied by the 
velocity field. 

In the light of these arguments, three possible pitfalls that can be encountered in the deter- 
mination of the local energy balance equation are described in the following. 

4.1 The isochoric process 

One could say that the convective derivative De/Di can be evaluated by assuming that the 
fluid element is undergoing an isochoric thermodynamic process [6]. Therefore the validity 
of Eq. ([9j) is extended, not only to the perfect gases, but to every fluid system. Then one 
would be led to Eq. @ instead of Eq. (|15|) . This reasoning is erroneous as the nature of the 
thermodynamic process undergone by the fluid element is implicitly defined by the equation 
of state assumed, i.e. Eq. ©. This equation of state is based on the hypothesis that the 
density may change only as a consequence of temperature variations. In other words, one is 
assuming that the set of stable equilibrium states available to the fluid element are, with a 
very good approximation, at constant pressure. 

4.2 The expansion-contraction work 

A possible misleading argument in the deduction of Eq. (|15|) is connected to the mechanical 
work term cr^ Dij. This term, as it is clearly explained in Chandrasekhar [6], must be simpli- 
fied according to the constraint satisfied by the velocity field, i.e. Eq. ([T]). On the other hand, 
in some textbooks (see, for instance, Kundu and Cohen p2]), a part of the mechanical work 
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term cr^Sy, namely the expansion-contraction work contribution, —pdui/dxi, is rewritten 
by forgetting Eq. ([T]) and by using the exact local mass balance instead, 

-»?r- £ & (20> 

OXi p Dt 

If one brings this term to left hand side of the local energy balance Eq. ((8j) , then one has to 
evaluate 

De p Dp 

P m-- P m (21) 

instead of pDe/Dt. Therefore, by following the thermodynamic argument described in Sec- 
tion [3l one would have 

p de — — dp = p c dT + p (3 dT = p ( c„ — - — | dT + p dT = pc p dT, (22) 
P \ P J 

where Eqs. (fTUI) . (fl~3j) and (fT4"j) have been used. On account of Eqs. (f20l) - (|22l) . one would be 

led to the formulation of the local energy balance expressed by Eq. ([7]) , instead of Eq. (I15p . 

However, this is a tricky procedure, since the velocity field evaluated through the Oberbeck- 

Boussinesq model is constrained to be solenoidal by Eq. (pQ). Then, even if it is based on an 

exact mass balance, Eq. (|20p cannot be coherently invoked within the Oberbeck-Boussinesq 

approximation . 

4.3 The pressure work 

One could evaluate the differential of the internal energy per unit mass, de, by using the 
definition of enthalpy per unit mass, h = e + p/p. Then, one has 

p 

pde = pdh — dp H — dp. (23) 
P 

If one assumes that the thermodynamic process undergone by the fluid element is isobaric, 
then one has 

dp = 0, dh = c p dT, dp = -p(3 dT. (24) 
Then, Eq. ([23D yields 

pde = pc p dT — p (3 dT = p ^c p — dT = pc dT, (25) 
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where Eq. has been used. Eq. (|25|) leads directly to Eq. (|15p . 

If one assumes that the thermodynamic process undergone by the fluid element is isochoric, 
then one has 

dp = 0, dh = c p dT + J dp. (26) 
The thermodynamic identity 

can be easily proved on the basis of the elementary thermodynamic differential relationships 
[14j . The complete proof can be found, for instance, in Ref. [13]. Therefore, by substituting 
Eqs. d2SD and ([27]) into Eq. ([23]), one obtains 

pde = pc p dT + (1 — (3T) dp — dp = pc p dT — f3T dp. (28) 

From Eq. (|28p . one justifies the relationship 



De DT nrTi Dp 



Obviously, Eq. ([29]) leads to Eq. (p. 

While the procedure based on Eqs. ([23]) - ([25]) and leading to Eq. (fT5|) is correct, the pro- 
cedure based on Eqs. ([231) . ([26]) - ([29]) and leading to Eq. ([6]) is completely unjustified for the 
following two reasons. 

• The process undergone by the fluid element cannot be considered as isochoric as already 
pointed out in Subsection 14. 11 

• If one has to suppose that the process is isochoric, the straightforward reasoning is to 
assume the general validity of Eq. Q, so that one is led to Eq. ([3]), as described in 
Subsection 14.11 

In other words, the misleading thermodynamic analysis behind Eq. ([6]) adds an unjustified 
assumption (an isochoric process) to a tricky procedure that replaces the simple conclusion 
De DT 

p Wt =pc "m (30) 



with the more complicated one expressed by Eq. (I29p . This circumstance induces some 
reflections on how the erroneous assumption of isochoric process can lead to ambiguous 
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results, depending on the procedure followed. On the contrary, the assumption of isobaric 
process leads exactly to the same result, i.e. Eq. (I15p . either if the reasoning stems from 
the evaluation of de or, as in the present subsection, if one works out the differential of the 
enthalpy per unit mass, dh. 

5 On porous media 

The remarks on the appropriate form of the local energy balance in the framework of the 
Oberbeck-Boussinesq approximation can be easily reformulated with reference to the theory 
of fluid saturated porous media. In fact, following Nield and Bejan [16], the volume-averaged 
energy equations for the solid and fluid phases can be written as 
dT 

(1 - tp) (Pcv)s^r = (!-¥>) k sV 2 T s + H (Tf — T s ) , (31) 
dTt dTf 

l p( pc ) f ^l + (pc) f V i ^- = <pk f V 2 T f + ^ + H(T s -T f ), (32) 

where V{ is the seepage velocity and the subscripts s and / denote the solid and the fluid 
phase, respectively. In Eq. (f32|) . ip is the porosity, while ^<3? is the viscous dissipation term 
obtained through a volume average of the term 2 p,Dij Dij appearing in Eq. (I15p . As is well 
known [17} 118] . the specific form of the term \x $ depends on the momentum balance model 
adopted for the fluid saturated porous medium. The inter-phase heat transfer coefficient 
H in Eqs. (|3ip and (|32p describes the thermal energy flow between the solid and the fluid 
phase. Since reference is made to the Oberbeck-Boussinesq approximation, all the solid and 
fluid properties in Eqs. (f3"Tj) and (|3"2"|) are evaluated at the reference temperature To. With 
respect to Nield and Bejan [16], Eq. (j32j) has been adapted on the basis of Eq. (fT5|) in order to 
include the thermodynamic coefficient c defined through Eq. (I13p . In cases of local thermal 
equilibrium between the solid phase and the fluid phase, T s = Tj = T, one can add Eqs. (f3T|) 
and (1321). so that one has 



dT dT 

(pc) m — + (pc) f Vi— = k m V 2 T + LL<5>, (33) 

where 

(pc) m = (l-<p) (pc v ) s + ip {pc) f , 

(34) 

k m = (1 - *p) k s + <pk f . 
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6 Conclusions 

The nature of the Oberbeck-Boussinesq approximation has been revisited in order to clar- 
ify some thermodynamic aspects connected to the formulation of the energy balance. The 
present analysis, motivated by the manifold formulation of this balance in the existing liter- 
ature, leads to the following conclusions: 

1. the energy balance for a fluid, within the Oberbeck-Boussinesq approximation, is given 
by the equation 

DT 

p c — = k V 2 T + 2 /x ®ij %, 
where c = c p — p/3 / p; 

2. the thermodynamic coefficient c coincides with the specific heat at constant volume, 
Cy, for a perfect gas; 

3. the thermodynamic coefficient c is definitely well approximated by the specific heat at 
constant pressure, c p , for a liquid; 

4. no pressure work term of the type f3 T Dp /Di must be introduced on the right hand 
side of the energy balance equation. 

Finally, the most appropriate form of the energy balance for buoyant flows in a fluid saturated 
porous medium has been discussed. 
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